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I. Introduction 

State feedback control system design is often related to reconstruction of the state vector by measurements of the output 
variable and the input signal of the open loop system. 

The reconstruction of the state vector is only possible by implication of state observer and the adaptive observation problem 
is related to observer synthesis with parameter estimator [6,7]. The matrices A and b or c (depending on the canonical form 
chosen for state space representation) are considered unknown. 

The parameters are being estimated and the unknown matrices are determined during the observation process and the state 
vector is reconstructed. 


The present paper investigates a non-recurrent algorithm for adaptive observation of single input single output (SISO) linear 
time invariant (LTI) discrete systems developed on the basis of the instrumental variables (IV) method [2], 

The parameters estimator built-in in the adaptive observer is based on a simplified calculation procedure which also includes 
inversion of the informative matrix [5], 


II. Problem Formulation 


The system investigated is presented in the state space with the following systems of equations: 


where: 


x(k + 1) = Ax(k) + bu(k), 
y(k) = c T x(k) + f(k) 
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The system order n is a-priori known, x(k) £ R n is the unknown current state vector, x(Q) £ R" is the unknown initial state 
vector, uQt) £ R l is a scalar input signal, y{k) £ R 1 is a scalar output signal, f(k) is an additive noise signal, a and b are 
unknown vector parameters. 


The following discrete transfer function corresponds to (1): 

yY(z) — h lz n ~ 1 +h 2 z n - 2 +"-+h n _ 1 z+h n 
z n —a n z n ~ 1 - a 2 z —a i 


(4) 
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The conformity between the elements b, of vector b in relation to the chosen phase canonical form for representation and the 
coefficients h,of the polynomial in the numerator of the discrete transfer function (4) are defined by the following 
expression[4]: 


Tb = h 


( 5 ) 


where: 
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The elements a, of vector a, in relation to the chosen phase coordinate canonical form for representation, are the coefficients 
of the polynomial in the denominator of the discrete transfer function (4) presented in reversed order and with an opposite 
sign. 

The problem to be solved is estimation of the unknown vector parameters a and b, the initial state vector x(0) and the 
current state vector x(k), k — 1,2 ,.... 

III. Solution 

3.1 Algorithm for adaptive observation based on the instrumental variables (IV) method 

The developed algorithm is based on the below calculation procedure: 

Step 7.Form the input-output data massive [3,8]: 

u x = [u(0) u(l) ••• u{N — 2)] 

yi = [y(o) y(i) ••• y(w-i)] 
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Where ^ll> ^21 , U 12 and U 22 are Toeplitz matricesand N — 3n + 21, l — 1,2,3,... 

Step 2.Calculate the sub-matrices G 11; G 12 , G 2 i, G 22 : 

Gu = + Y 21 Y 21 ; G 12 = K[,G 12 + Y T 2l U 22 , G 21 = G[ 2 K u + C 22 = G[ 2 G 12 + Gl 2 G 22 

Step 3. Calculate the covariance matrix C: 

M l + M 1 G 12 M 2 G 21 M 1 : —M 1 G 12 M 2 

—M 2 G 21 Mi : M 2 

where: 

Mi — G 11 1 ;M 2 = (G 22 — G 21 M 1 G 12 ) 1 

Step 4. Calculate the vectors h and a by using the following vector-matrix system and form the estimated system matrix A: 

\Yhy 2 +Y T 21 y 3 ] 


c = 


v — c 


h — [h 1 h 2 


L U T 12 y 2 + U T 22 y 3 \ 

h n ¥ = \Pn+i Pn +2 "• p 2 „] 7 ';a=[a 1 a 2 ••• a n ] T = [~p n -f n _i 

o : /„_r l 

a T 


-PiY 


A = 


Step 5. Calculate vector b estimation by implementing the following linear algebraic system of equations: 

Tb — h 


is a lower triangular Toeplitz matrix. 


Step 6. Estimate the initial state vector x 0 : 

*o = (D T D)~ 1 D T (y 1 - QuJ = [x Q1 x 02 ••• x 0n 

(applicable only in case that det(D T D) ^ 0), where: 
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Step 7. Calculate the estimate of the output variable y(fc): 

x(k + 1) = Ax(k) + bu(k), 
y(k ) = c T x(k ) 


*(o) = *o 

k = 0,1,2,... ,1V — 1 


F — A — gc T 


Step tf.For the instrumental matrices E u , V 21 : 
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Step ^.Recalculate the submatrices G n and C 12 : 

Gn = V t u Y u + ^y 21 ; G 12 = V t u U 12 + G 22 

Step 10. Recalculate the parameters vector p: 

Mi — G 11 1 ;M 2 = (G 22 — G 21 M 1 G 12 ) 1 
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Step 77.Repeat steps 7 to 10 four times 

Step 72.Estimate the current state vector x(k): 

x(k + 1) = Fx(k) + ftu(fe) + gy(k). 5c(0) = 5c 0 ; F = A — gc T 

Vector g can be easily obtained by solving the pole placement problem (PPP) also known as pole assignment problem 
(PAP). The synthesis of vector g must take into consideration the following options: the eigen values of the matrix F should 
be zeros or should be spread into the unit circle closer to the origin than the matrix A eigenvalues. The implementation of the 
above listed options ensures good dynamic characteristics of the observer synthesized. 

IV. Simulation results 


The simulation is held in MATLAB programming environment under the below conditions: 

> For a given transfer function of the system under investigation, with input signal u(k) and the respective output 
signal y(k); 

> Added colored noise signal f(k) is applied to the system output; 

> The input signal u(k) and the noise-corrupted output signal y(/c) are used as an input data for the observation algorithm; 

> Based on the input-output data massive, the algorithm developed calculates the estimates of the open loop system and the 
state vector. 


The discrete transfer function of the system investigated, used for the simulations is presented as follows: 


W (z) = 


0.6z _1 + 0.56z -2 + 0.2125z -3 + 0.308z“ 4 + 0.5488z“ 5 + 0.7221z" 6 
1 - 1.4z _1 + 0.7875z" 2 - 0.2275z- 3 + 0.035525z- 4 - 0.002835z- 5 + 0.00009z- 6 


and its corresponding representation in state space: 
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The matrix A eigenvalues are obtained by using eig(.) function in MATLAB: 

eig{A) = [0.4 0.3 0.25 0.2 0.15 0.l] r 

Pseudo-random binary sequence (PRBS) is used as an input signal u(k) which is generated in MATLAB by using the 
following functions: u — [sign{randn{ 127,1)) j * 10. 


The output signal y(/c) is noise-corrupted by adding a color noise /(/c). The colored noise is obtained by filtering of white 
noise through filter with the following transfer function: 


W f (z) = 


1 

1 - 1.4z" 1 + 0.7875z- 2 - 0.2275z“ 3 + 0.035525z- 4 - 0.002835z- 5 + 0.00009Z" 6 


The noise level g is calculated mathematically by division of the noise standard deviation <je to the output signal standard 
deviation a y in accordance with the equation given below: 

r\ = XlOO = (0 h- 10)% 

°y 

Vector a estimation error e a , vector b estimation error e b and the state vector x(k) estimation error e x are relative mean 
squared errors (RMSE) and could be determined by the following equations: 


e a (k) = - 
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e b (fc) = - 
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The results for the case of noise-free output signal (i.e. f(k) — 0) and l — 0 (i.e. minimum number of the input-output 
measurements: N = 3n — 18) are shown in Fig.l. In accordance with the initial settings the algorithm will start working at 
the 18 th step of the calculations and in this case in particular the observation errors e a (k), e b (J<) and e x (k) are zeros. 



k 

Figure 1. RMSE for the case of noise-free output signal 
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In the case of noise-corrupted output signal an experiment is carried out for noise level 77 = 10.018% and l — 40 (i.e. 
N — 3n + 21 — 98). The results are shown in Fig.2. Under the above described initial settings the algorithm will start 
working at the 98 th step of the calculations and the RMSE are as follows: e a (fc) < 0.033, e b (fc) < 0.010, e x (fc) < 0.065. 
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Figure 2. RMSE for the case of noise-corrupted output signal, tj = 10.018%, N = 3n+ 21 = 98 

An additional experiment is carried out for the provisions of noise-corrupted output signal analysis for noise level r/ — 
10.014% and l — 100 (i.e. N — 3n + 21 = 218). The results show that the algorithm starts working at the 218 th step of the 
calculations and the RMSE are as follows: e a (fc) < 0.017, e b (fc) < 0.0057,e x (fc) < 0.024, under the condition:218 < k > 
400 => e x {k) = 0.01 (see Fig.3). 
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FIGURE 3. RMSE FOR THE CASE OF NOISE-CORRUPTED OUTPUT SIGNAL,^ = 10. 014 %,N = 3n + 21 = 218 
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The simulation results delivered and graph analysis show that with the increase in number of the input-output measurements 
(AO, the algorithm invariance to added noises increases proportionally however the time needed for collection of initial 
information rises. 


V. Conclusions 

The algorithm suggested for open loop system parameter estimations which serve as a basis for further reconstruction of the 
current state vector implements the method of the instrumental variables excluding the zero iteration which only uses the 
least squares method (steps 1 to 4 of the suggested calculation procedure). 

The algorithm proposed estimates as well the initial state vector x 0 which allows the forming of the instrumental variables 
matrix even for nonzero initial conditions. 

The results delivered show that the number of the input output data measurements (AO is of high significance in relation to 
accuracy of estimations in the case of noise corrupted output. The highest accuracy is to be expected for highest counts of N 
(see Fig.2 and Fig.3). 

The method of the IV method gives best results in case of estimation of a-priori collection of data [3,8], however in relation 
to the closed loop system the added noise f(k) is transferred to the input signal through the feedback channel. Thus 
invariance between the instrumental matrices and the added noise is not possible; it is only possible that the estimates are 
unbiased and significant in presence of a white noise however the real systems do not allow such solution of the problem. 
The used of IV method for investigation of the closed loop system is only applicable if additional input signal is implemented 
[8], For this reason the implementation of the algorithm suggested above is not recommended for closed systems 
implications. 

The algorithm for adaptive observation based on the instrumental variables (IV) method introduced in the present paper is 
developed on the basis of non-recurrent method which ensures the convergence of the iterative procedure [1,2], 

The most positive feature of this algorithm however is related to the method used for informative matrix formation. It is 
formed through the four sub-matrices ^ll< ^21 , U 12 and U 22 which reduces the calculation complexity of the procedure for 
inversion of matrix G formed by the sub-matrices G n , G 12 , G 2 \, G 22 . Independently of the N count in numbers for the 
estimation of the coefficients h t anda,is only needed to be inverted the matrices G n and ( G 22 — G 21 M 1 G 12 ) which are 
always guaranteed (n X n) dimensions. In all other cases this procedure is related to inversion of a matrix at least (N — n) X 
(N — n) dimensional. 
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